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Abstract: Two novel classes of spherical invisibility cloaks based on 
nonlinear transformation have been studied. The cloaking characteristics are 
presented by segmenting the nonlinear transformation based spherical cloak 
into concentric isotropic homogeneous coatings. Detailed investigations of 
the optimal discretization (e.g., thickness control of each layer, nonlinear 
factor, etc.) are presented for both linear and nonlinear spherical cloaks and 
their effects on invisibility performance are also discussed. The cloaking 
properties and our choice of optimal segmentation are verified by the 
numerical simulation of not only near-field electric-field distribution but 
also the far-field radar cross section (RCS). 
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1. Introduction 

The use of coordinate transformation to control electromagnetic fields yj has been receiving 
extensive attention ||2l[3]|4l. This approach was generalized from the cloaking of thermal con- 
ductivity [f5|, and was further widely applied in electromagnetics and acoustics |[l][6l, which 
provides new ways to conceal passive/active objects Q [S) within their interiors invisible to 
external illuminations. The fundamental idea is that Maxwell's equations are invariant under a 
coordinate transformation if the material properties (electric permittivity and magnetic perme- 
ability) are altered appropriately; i.e., a specific spatial compression is equivalent to a variation 
of the material parameters in flat (original) space. Cloaks in regular shapes (e.g., spherical HI, 
square lH, cylindrical fTOl , or elliptic ||TT|) have been proposed based on the construction of 
explicit transformation matrices. There are many works devoted to the investigation of 2D cylin- 
drical/elliptic cloaks approximated by, for example, simplified material parameters or equations 
lfT2l[T3][T4l[T5l . or designed for limited incident angles lfT6l . Experimental investigations of 2D 



cloaks have demonstrated significant reductions in the scattering cross section, albeit for narrow 
bandwidths in the microwave regime and for objects at most a few wavelengths in diameter so 
far iflTl . Broader bandwidths have been demonstrated for cloaking an object on a ground plane, 
but again only for a few wavelengths diameter [18 |. To design irregularly shaped cloaks, both 
analytical ||T9| and numerical methods ll20l |2TI have been proposed. In what follows, we will 
confine our discussion within the area of 3D electromagnetic spherical cloaks. 

Inspired by the classic spherical cloak Ul based on a linear coordinate transformation, the ex- 
pressions of electromagnetic fields were explicitly presented in terms of spherical Bessel func- 
tions along the lines of Mie scattering theory |f7|, and it again was demonstrated that the external 
incident wave cannot interact with the cloaked object in a perfect spherical cloak. Linear coordi- 
nate transformations correspond to a transverse/radial permittivity ratio {anisotropy ratio 1221 ) 
et/er = I {r^Ri where R\ is the radius of the cloaked region. The subtle point of the singu- 
larity in the coordinate transformation at the innner surface of the cloak was analytically shown 
to correpond to surface voltages |8| in spherical cloaks, leading to an explicit physical expla- 
nation for why the wave cannot leave the cloaked region and external waves cannot enter that 
region. These works on spherical cloaks considered linear coordinate transformations, while 
more recently nonlinear ("high-order") transformations have been considered in order to obtain 
more degrees of freedom in designing the material parameters ||23ll24l . In this paper, we con- 
sider the utility of nonlinear coordinate transformations to improve the performance of the cloak 
when it is approximated (segmented) into a sequence of piecewise-homogeneous layers. If the 
anisotropy ratio differs from r^/(r — it will lead to quite complicated formulations for 
the field expressions, which cannot be treated by the previous method for a position-dependent 
anisotropy ratio |7:| or by the method for a constant anisotropy ratio 1251 |26l . Here, in order to 
study the proposed nonlinear-transformation spherical cloaks, each anisotropic cloaking shell 
is discretized into multiple alternating homogeneous isotropic coatings, not only for the ease of 
fabrication but also for the computation convenience in dealing with arbitrary anisotropy ratios. 

In general, practical implemenations of passive invisibility cloaks by coordinate transfor- 
mation are limited by several factors, in addition to the bandwidth limitations (from material 
dispersion) mentioned above. First, there are fabrication imperfections and the finite size of the 
subwavelength components of the metamaterials. Second, there is material absorption, which 
is especially challenging when an extreme resonant response is required of the material in or- 
der to obtain very large or very small permittivity/permeability. Third, there is the difficulty of 
fabricating continuously varying anisotropic materials. In this paper, we focus only on the third 
issue. A common approximation for continuously varying materials is to segment them into 
piecewise-homogeneous layers, and we show that the choice of coordination transformation 
has a significant impact upon the success of this approximation. We argue that the segmenta- 
tion is best performed in virtual space rather than in the physical (transformed) coordinates. 
We also demonstrate that the anisotropic materials can be further approximated by a sequence 
of isotropic -material layers. The optimal selection of material parameters, nonlinear transfor- 
mation, and segmentation (in virtual space), in order to achieve a low scattering cross sectikon 
over a wide range of observation angles, is found for nonlinear spherical cloaks and verified by 
numerical results. 

2. Preliminaries 

The configuration of the spherical cloaking structure is shown in Fig. [T| The inner region (r < 
^i) is a perfect electric conductor (PEC) and the intermediate region {Ry < r < Rt) is filled 
by the nonlinearly transformed spherical cloak, characterized by e(r) and jil(r), which will be 
discussed below. The electric field is polarized along the x axis and propagating along the z axis. 
By applying the decomposition method and separation of variables ||22L the equation for the 
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Fig. 1 . The geometry of a spherical cloaking structure. 



radial component of the field potentials becomes 

n(n + 1) 
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where the subscripts e and m denote the electric and magnetic anisotropic ratio, respectively, 
and kt ~ (O^/J^Ft EU- For £ and fl of perfect linear spherical cloaks |[T|, i.e., Ae = A,„ = 
/{r — Rif', Eq. ([T]) is reduced to 
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Thus, the radial component f{r) can be solved in a way similar to that for isotropic materials, 
except for the change in the argument of resultant Bessel/Hankel functions. However, given a 
set of e and /i derived from a certain transform, the anisotropy ratio may not be r^/ {r — RiY 
anymore, and then the radial component cannot be solved explicitly in the same way. In this 
situation, we could approximate the original inhomogeneous anisotropic cloak by many thin, 
homogeneous anisotropic coatings, and the diffraction problem can thus be solved in terms of 
analytical Bessel/Hankel functions satisfying boundary conditions at each interface. Neverthe- 
less, the requirement of anisotropic materials still remains. Alternatively, to further alleviate the 
restriction in material complexity, the original inhomogeneous anisotropic cloaking materials 
can be approximated by the limit of many thin, concentric, homogeneous, isotropic coatings, 
forming an effective anisotropic medium, which will be discussed in this section. 

First, we consider a general spherical cloak characterized by £ = e,.(r)rr + e,(r)(00 + p^) 
and 11 = jir{r)7?+ jj.,{r){99 + 00), in which er{r) = IJ-ri'') and £,(r) = ji,{r) are position- 
dependent in general. It is then divided into M initial-layers (anisotropic) but the thickness of 
individual initial layers may or may not be identical, depending on the transformation. Given a 
coordinate transformation function r' = f{r) between the virtual space (i.e., ii'(r'), < r' < R2) 
and the compressed space (i.e., i2(r), Ri < r < R2), the stepwise segmentation in physical space 
(r) is desired to mimic the transformation function as well as possible and we also desire that 
the segmentation not be too complicated. We find that equally dividing the virtual space (r') 
into M initial-layers will make the segmentation on the physical space (r) "self-adaptive" in a 
simple way — it automatically uses a finer discretization in regions of physical space where the 
anisotropic materials are varying more rapidly — which will result in better invisibility perfor- 



(a) anisotropic and inhomogeneous 



(b) isotropic and homogeneous 



Fig. 2. The discretization of a general anisotropic splierical cloak (a) into an equivalent 
isotropic coated structure (b). An illustrative example for the conversion from the n-th 
anisotropic initial-layer (width is r„ — r„_i) into its two sub-layers (isotropic) with equal 
thickness {r„ — r„_ i ) /2 has been shown. 

mance. This conversion is illustrated in Fig.|2l Throughout this paper, we apply the same con- 
dition, i.e., the segmented layers in r' are of equal thickness, which represents a good choice of 
segmentation in the initial-layers in the cloak shell Ri < r < Ri- By projecting the segmentation 
in r', = R2 ■ n/M onto the physical r, one has r„ = f^^{f^„), (n=l, 2,..., M). Thus, the geometry 
of every initial- layer in Fig. 12 a) is determined. 

Subsequently, we mimic each spherically anisotropic initial layer by a pair of effective 
isotropic sub-layers (type-A and type-B) with equal thickness at (r„ — r„_i)/2, as shown in 
Fig. |2b). The material parameters of type-A and type-B isotropic dielectrics can be inferred 
and derived from the result for radial conductivity by Sten |f27|; 



Now, the original spherical cloak turns to be a concentric isotropic coatings. Throughout this 
paper, M = 40 and R2 = 2Ri = 2A are chosen in order that the validity of Eq. (O can be 
maintained, i.e., sub-layers have to be sufficiently thin. Alternatively one could also choose a 
smaller M associated with the corresponding parameters determined by optimization methods 
and then subdivide each anisotropic layer into many more than two isotropic sub-layers. 

When converting the «-th initial layer into its pair of effective isotropic sub-layers, we need 
to pick a specific radial position for those e, and e, on the right-hand side of Eq. Q within 
the initial layer in order to determine corresponding parameters of the two isotropic dielectrics 
(Ea, Ha) and (eg, ^b) on the left-hand side of Eq. ^ and in Fig. I^b). According to 1291 . 
such a discretization mechanism using r ~ r„ in Eq. ^ for the «-th initial layer will give a 
good compromise for both forward and backward scatterings, while retaining good invisibility 
performance. 

3. Effects of Nonlinear Transformation in Nonlinear Spherical Cloaks 

Pendry's spherical cloak lH] is a Unear one, and hence it is a straight line if one plots the trans- 
formation function r' against the physical radius r. In what follows, we introduce two classes 
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of nonlinear-transformation spherical cloaks, and we discuss how to restore and improve the 
invisibility performance after discretization by choosing a proper nonlinear transformation and 
by choosing a suitable compensation scheme while discretizing the original spherical cloak into 
multilayer isotropic structures. The two types of spherical cloaks considered here are classified 
in terms of the negative (i.e., concave-down) or positive (i.e., concave-up) sign of the second 
derivative of the transformation function. It is important to reiterate that all three designs — 
linear, concave-up, and concave-down — are perfect cloaks for the exact inhomogeous design, 
and we are only considering the breakdown of invisibility when the design is discretized into 
homogeneous layers. 

3.1. Concave-Down Nonlinear Transformation 

To compress a sphere of air at the radius R2 in Q.' (original) space into a shell at the region 
Ri < r < R2 in Q. (compressed) space, we propose a class of prescribed functions 
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where x denotes the degree of the nonlinearity in the transformation. When x is very small in 
Eq. (HI, the curves are difficult to distinguish and all approach to the same limiting case when 
x^O: 

RoLnlr/Ri] 
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Thus the parameters e, Jl in the transformed coordinates can be written in term of e', Jl' in 



the original space by 

e(r) =Ae'{r')A'^/det{A), J[{r) ^ AJl' {r')A^ /det{A) 
where A is the Jacobian matrix with elements defined as A^^, = dr/^/ dr'. 
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Fig. 3. The concave-down nonlinear transformation. When x < 0.1, all mapping curves 
are overlapping, meaning they effectively lead to the same performance. When x becomes 
large, the curve will nearly become a step function over IX < r <2X. 



Note that all curves belonging to those transformation functions in Eq. (HI have negative 
second derivative with respect to the physical space r, and we term this class of transformations 



the concave-down nonlinear transformation. The nonlinear transformation function in Eq. (|4]i 
only depends on the radial component r in the spherical coordinate system (r, 9, (p). Thus it is 
easy to find that the Jacobian matrix A is diagonal. Considering that the original space is filled 
with air (e' —Jl' — 1), Eq. (|6]l can be rewritten as 

e = /r = dmg[Xl Xl, Xl]/XrXeX^ = diag[^, ^] (7) 

AgAA A,-A|f, A,-A0 



where 



^r=^ = — ^TT (8a) 



:A. = ^= ^ (8b) 



denoting three principal stretches of the Jacobian matrix. Finally, the desired parameters in the 
compressed space (Ri < r < Rj) are shown to be 

R^+Ur-^-R^y 



xR] {R^2 - ^1 ; 



xR]Rl+' 



ee=^e = s,=^, = j^^-^^. (9b) 

Given such radial and transversal parameters, it is obvious that Eq. ^ cannot be simply solved 
due to the complicated anisotropy ratio Ae,,„ ~ et/£r= y?'R?^ j (r^ — )^ in our current case. In 
addition, one can also find that when the virtual space is equally discretized using a concave- 
down nonlinear transformation, the segmentation of the initial layers in Fig.|2]in the physical 
space r is 



r„(x) 



R\-R\ n 



R\ M 



■Rx, n=l, 2,...,M. (10) 



Here, we will study the effects of the nonlinear factor x in the near-field and far-field of the 
discretized nonlinear-transformation spherical cloaks. We fix M=40, and consider the cases of 
x=0.1, x=l, x=4, and x=10 in concave-down nonlinear cloaks, shown in Fig.|4] It can be seen 
from Figs. |4|a-d) that when x increases, the magnitude of electric field increases significantly 
inside the cloak. This is because more energy is guided towards the inner boundary of the 
PEC core, which in turn makes the cloaked PEC more visible to external incidences. To prove 
that the large electric fields only occur in the region Ri < r < R2 and to study the effect of 
X on individual near-field perturbations more explicitly, only the fields outside the cloak are 
presented in Figs. Ue-f). One can again confirm that, for discretized concave-down nonlinear 
spherical cloaks, smaller x leads to smaller disturbance in electric fields in the outer space. In 
particular, the degradation in the invisibility in the forward direction is proportional to the value 
of X (see Fig.Htg) and Fig.Sfh)). 

In Fig. |5] we consider the far fields of those cases of Fig. |4] For comparison purposes, we 
introduce two more curves, PEC and Pendry's cloak. It should be noted that the curve corre- 
sponding to Pendry's cloak is a straight line, and the same segmentation procedure as discussed 
above is applied to it, which is the same as equally dividing it into M initial-layers in the physi- 
cal space (see Fig. O. From Fig. |5] it can be seen that for concave-down nonlinear transforma- 
tions, the cloaking property is better retained when x is small. Except for the forward direction. 
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Fig. 4. Near-field interaction in the presence of the proposed nonlinear spherical cloak with 
different values of x for the concave-down class. The inner region < r < Ri is filled by 
PEC. M is set to be 40. Frequency is 2GHz. The total electric fields are plotted only in the 
region r > Ro in (e-h). 
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Fig. 5. Bistatic RCS of the concave-down nonlinear cloaks. M = 40. 



concave-down nonlinear transformation with small x can achieve even lower RCS over a wide 
range of observation angles. As x keeps increasing, their RCSs are increasing dramatically and 
can be larger than that of a uncloaked PEC core (e.g., x = 10). No matter how large x becomes, 
the suppressed backward scattering is still maintained |71 though it is not exactly zero due to the 
discretization, but the forward scattering will become higher and higher, which is not desired 
in the cloaking application. Those far-field characteristics are consistent with the perturbation 
in corresponding near-field patterns in Fig.|4] 

3.2. Concave-Up Nonlinear Transformation 

We can also propose another class of prescribed functions 

= (11) 

It is found that when x ^ 0, it also approaches to the same limit as Eq. (|5]l. In the same manner, 
we can obtain the desired parameters in the compressed space {R\ < r < R2) 
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Note that the case x = 1 is exactly Pendry's linear spherical cloak. 

It can be seen in Fig.|6]that all curves belonging to Eq. (fTTl i have positive second derivatives, 
and we term this class of transformations as the concave-up nonlinear transformations. Also, it 
is found that the hmiting case of Eq. (|5]l, i.e, x 0, is the dividing line between the convace- 
down and concave-up classes. Since x <C lin Eq. (fTTT i is actually very close to the situation of 
X ^ 1 for concave-down nonlinear transformations, we will not revisit that case here. Instead, 
we only consider x = 1, x = 4, and x = 10 concave-up transformations [Eq. (fT2li1. and the 
segmentation for concave-up nonlinear spherical cloaks in physical space r turns out to be 
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We are interested to see which class will give rise to better invisibility performance, after 
discretization, if all the other conditions are the same. We still plot the near field of concave- 
up cloaks at "x=l", "x=4", and "x=10". Comparing Fig. |4] with Fig. |7] one can see that: (1) 
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Fig. 6. The concave-up nonlinear transformation. Wiien x = 1, it is exactly Pendry's linear 
cloak, which can be verified by inserting x = I into Eq. l ll2t . When x gets extremely large, 
the mapping curve will nearly become a sharp impulse when r^2X. 



the peak amplitude of the electric field is proportional to x in both classes; (2) for each value 
of X, concave-up transformation cloaks have lower peak amplitudes than concave-down trans- 
formation cloaks, i.e., the perturbation inside the cloak region is smaller; (3) the perturbation 
outside the cloak of concave-up transformation cloaks is also lower than that of concave-down 
transformation cloaks in all cases of x; (4) the invisibility performance is better maintained for 
concave-up transformations even when x^ 10 or even larger. Hence, spherical cloaks based on 
concave-up nonlinear transformations exhibit better invisibility after discretization. 

Fig. [U exhibits some interesting features which distinguish it from the far fields of concave- 
down transformations in Fig. |5] Over a wide range of angles, the invisibility performance of 
spherical cloaks based on concave-up nonlinear transformations are more stable under varia- 
tions in X. More importantly, when x increases (e.g., x ~4 or x ~ 10), though the RCS will be 
a bit larger than Pendry's linear cloak at most angles, the RCS near the forward direction will 
be reduced greatly compared with that of classic linear cloak. Such far-field phenomena are 
also connected with the near-field patterns. It is important to point out that the increase of x in 
concave-up transformations [see Figs.|7];a-c)] will push the "hot" areas (where the electric field 
is very high) further and further away from the spherical PEC core, so the induced shadow and 
the far-field pattern become stable with increasing x. However, the increase of x in the concave- 
down transformations has little effect in shifting away the "hot" positions from the core [see 
Fig. Sta-d)], which results in larger interaction with the PEC core. Hence, the RCS reduction 
in the far field is degraded with the increase of x in the concave-down class as shown in Fig. |5] 

4. Justification of Improved Segmentation 

It has been pointed out that usually we require a large value of M during the discretization pro- 
cess, otherwise the conversion scheme in Fig.|2]will lose its accuracy. We also find that a nonlin- 
ear transformation with a given x of its concave-up class will yield a better invisibility by taking 
the segmentation as Eq. (T3[ in the physical space (this is also the only way to determine sub- 
layers of type-A and type-B dielectrics in practice). Here the justification is given. The concave- 



up class with x = 4, which is characterized by the mapping curve r (4) = / h {r/Ri) —I 

is considered as an example. 

In Table I, different segmentation r„ (x) sets (see Eq. [TJt for real radius r are applied to 
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Fig. 7. Re[EjQj^j] on x-z plane for the proposed concave-up nonlinear transformation cloaks 
at X = 1, X = 4 and x = 10, respectively. Figs. |7{a-c) present total electric fields in all 
areas while Figs.|7te-f) show the total fields only outside the cloak (r > R2). M = 40 and 
/ = IGHz- In order to show the disturbance in surrounding outer space, the plot range is 
larger (from -0.6m to 0.6m) in the right column. 
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Fig. 8. Bistatic RCS of the concave-up nonlinear cloaks. 



Table 1. The total scattering cross section (unit: dBsm) for concave-up nonlinearly trans- 
formed spherical cloaks (x = 4) under various segmentations in physical space. 



r„{x) in Eq. (13) 


'•«(!) 


rni2) 


rni4) 


rni6) 


r„(10) 


r'(4) in Eq. (11) 


-3.3 


-3.5 


-18.9 


-3.8 


-3.1 



. * Note that r„{l) for concave-up cloaks corresponds to the case of equally dividing the physical space. 



r'(4), e.g., the sets of jc ^ 1 (r„(l)), jc = 2 (r„(2)), jc = 6 (r„(6)) andx 10 (r„(10)). The total 
cross sections corresponding to those four sets are compared with the best segmentation set of 
its own r„(4), i.e., equally dividing the virtual space of jf = 4 into 40 initial-layers and then 
projecting onto physical space. We find that, for concave-up class x = 4, r„(4) is indeed its 
optimal segmentation among these possibihties. 

5. Conclusion 

In this paper, we discussed two novel classes of nonlinear-transformation-based spherical 
cloaks. An approximation mechanism for such classes of spherical cloaks by concentric 
isotropic sub-layers is proposed. The optimal segmentation in the virtual space is presented, 
and the role of nonlinearity in the coordinate transformation is discussed in order to provide 
better invisibility performance for the proposed classes of nonlinear spherical cloaks in both 
the near and far field. We find that concave-up nonlinear-transformation spherical cloaks exhibit 
better and more stable invisibility performance than the corresponding segmentation of linear- 
transformation or concave-down cloaks. An interesting topic for future exploration would be to 
consider a wider class of segmentation schemes, for example by using the segmentation pre- 
sented here as the starting point for a numerical optimization (e.g., Ref. Il28l for cylindrical 
cloaks) of variable segment thicknesses and material parameters. 
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1. Introduction 

The use of coordinate transformation to control electromagnetic fields yj has been receiving 
extensive attention ||2l[3]|4l. This approach was generalized from the cloaking of thermal con- 
ductivity [f5|, and was further widely applied in electromagnetics and acoustics |[l][6l, which 
provides new ways to conceal passive/active objects Q [S) within their interiors invisible to 
external illuminations. The fundamental idea is that Maxwell's equations are invariant under a 
coordinate transformation if the material properties (electric permittivity and magnetic perme- 
ability) are altered appropriately; i.e., a specific spatial compression is equivalent to a variation 
of the material parameters in flat (original) space. Cloaks in regular shapes (e.g., spherical HI, 
square lH, cylindrical fTOl , or elliptic ||TT|) have been proposed based on the construction of 
explicit transformation matrices. There are many works devoted to the investigation of 2D cylin- 
drical/elliptic cloaks approximated by, for example, simplified material parameters or equations 
lfT2l[T3][T4l[T5l . or designed for limited incident angles lfT6l . Experimental investigations of 2D 



cloaks have demonstrated significant reductions in the scattering cross section, albeit for narrow 
bandwidths in the microwave regime and for objects at most a few wavelengths in diameter so 
far iflTl . Broader bandwidths have been demonstrated for cloaking an object on a ground plane, 
but again only for a few wavelengths diameter [18 |. To design irregularly shaped cloaks, both 
analytical ||T9| and numerical methods ll20l |2TI have been proposed. In what follows, we will 
confine our discussion within the area of 3D electromagnetic spherical cloaks. 

Inspired by the classic spherical cloak Ul based on a linear coordinate transformation, the ex- 
pressions of electromagnetic fields were explicitly presented in terms of spherical Bessel func- 
tions along the lines of Mie scattering theory |f7|, and it again was demonstrated that the external 
incident wave cannot interact with the cloaked object in a perfect spherical cloak. Linear coordi- 
nate transformations correspond to a transverse/radial permittivity ratio {anisotropy ratio 1221 ) 
et/er = I {r^Ri where R\ is the radius of the cloaked region. The subtle point of the singu- 
larity in the coordinate transformation at the innner surface of the cloak was analytically shown 
to correpond to surface voltages |8| in spherical cloaks, leading to an explicit physical expla- 
nation for why the wave cannot leave the cloaked region and external waves cannot enter that 
region. These works on spherical cloaks considered linear coordinate transformations, while 
more recently nonlinear ("high-order") transformations have been considered in order to obtain 
more degrees of freedom in designing the material parameters ||23ll24l . In this paper, we con- 
sider the utility of nonlinear coordinate transformations to improve the performance of the cloak 
when it is approximated (segmented) into a sequence of piecewise-homogeneous layers. If the 
anisotropy ratio differs from r^/(r — it will lead to quite complicated formulations for 
the field expressions, which cannot be treated by the previous method for a position-dependent 
anisotropy ratio |7:| or by the method for a constant anisotropy ratio 1251 |26l . Here, in order to 
study the proposed nonlinear-transformation spherical cloaks, each anisotropic cloaking shell 
is discretized into multiple alternating homogeneous isotropic coatings, not only for the ease of 
fabrication but also for the computation convenience in dealing with arbitrary anisotropy ratios. 

In general, practical implemenations of passive invisibility cloaks by coordinate transfor- 
mation are limited by several factors, in addition to the bandwidth limitations (from material 
dispersion) mentioned above. First, there are fabrication imperfections and the finite size of the 
subwavelength components of the metamaterials. Second, there is material absorption, which 
is especially challenging when an extreme resonant response is required of the material in or- 
der to obtain very large or very small permittivity/permeability. Third, there is the difficulty of 
fabricating continuously varying anisotropic materials. In this paper, we focus only on the third 
issue. A common approximation for continuously varying materials is to segment them into 
piecewise-homogeneous layers, and we show that the choice of coordination transformation 
has a significant impact upon the success of this approximation. We argue that the segmenta- 
tion is best performed in virtual space rather than in the physical (transformed) coordinates. 
We also demonstrate that the anisotropic materials can be further approximated by a sequence 
of isotropic -material layers. The optimal selection of material parameters, nonlinear transfor- 
mation, and segmentation (in virtual space), in order to achieve a low scattering cross sectikon 
over a wide range of observation angles, is found for nonlinear spherical cloaks and verified by 
numerical results. 

2. Preliminaries 

The configuration of the spherical cloaking structure is shown in Fig. [T| The inner region (r < 
^i) is a perfect electric conductor (PEC) and the intermediate region {Ry < r < Rt) is filled 
by the nonlinearly transformed spherical cloak, characterized by e(r) and jil(r), which will be 
discussed below. The electric field is polarized along the x axis and propagating along the z axis. 
By applying the decomposition method and separation of variables ||22L the equation for the 
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Fig. 1 . The geometry of a spherical cloaking structure. 



radial component of the field potentials becomes 

n(n + 1) 



— A 



}fM = o, 



(1) 



where the subscripts e and m denote the electric and magnetic anisotropic ratio, respectively, 
and kt ~ (O^/J^Ft EU- For £ and fl of perfect linear spherical cloaks |[T|, i.e., Ae = A,„ = 
/{r — Rif', Eq. ([T]) is reduced to 



«(« + 1) 



{r-Ri 



}f(r)^0. 



(2) 



Thus, the radial component f{r) can be solved in a way similar to that for isotropic materials, 
except for the change in the argument of resultant Bessel/Hankel functions. However, given a 
set of e and /i derived from a certain transform, the anisotropy ratio may not be r^/ {r — RiY 
anymore, and then the radial component cannot be solved explicitly in the same way. In this 
situation, we could approximate the original inhomogeneous anisotropic cloak by many thin, 
homogeneous anisotropic coatings, and the diffraction problem can thus be solved in terms of 
analytical Bessel/Hankel functions satisfying boundary conditions at each interface. Neverthe- 
less, the requirement of anisotropic materials still remains. Alternatively, to further alleviate the 
restriction in material complexity, the original inhomogeneous anisotropic cloaking materials 
can be approximated by the limit of many thin, concentric, homogeneous, isotropic coatings, 
forming an effective anisotropic medium, which will be discussed in this section. 

First, we consider a general spherical cloak characterized by £ = e,.(r)rr + e,(r)(00 + p^) 
and 11 = jir{r)7?+ jj.,{r){99 + 00), in which er{r) = IJ-ri'') and £,(r) = ji,{r) are position- 
dependent in general. It is then divided into M initial-layers (anisotropic) but the thickness of 
individual initial layers may or may not be identical, depending on the transformation. Given a 
coordinate transformation function r' = f{r) between the virtual space (i.e., ii'(r'), < r' < R2) 
and the compressed space (i.e., i2(r), Ri < r < R2), the stepwise segmentation in physical space 
(r) is desired to mimic the transformation function as well as possible and we also desire that 
the segmentation not be too complicated. We find that equally dividing the virtual space (r') 
into M initial-layers will make the segmentation on the physical space (r) "self-adaptive" in a 
simple way — it automatically uses a finer discretization in regions of physical space where the 
anisotropic materials are varying more rapidly — which will result in better invisibility perfor- 



(a) anisotropic and inhomogeneous 



(b) isotropic and homogeneous 



Fig. 2. The discretization of a general anisotropic splierical cloak (a) into an equivalent 
isotropic coated structure (b). An illustrative example for the conversion from the n-th 
anisotropic initial-layer (width is r„ — r„_i) into its two sub-layers (isotropic) with equal 
thickness {r„ — r„_ i ) /2 has been shown. 

mance. This conversion is illustrated in Fig.|2l Throughout this paper, we apply the same con- 
dition, i.e., the segmented layers in r' are of equal thickness, which represents a good choice of 
segmentation in the initial-layers in the cloak shell Ri < r < Ri- By projecting the segmentation 
in r', = R2 ■ n/M onto the physical r, one has r„ = f^^{f^„), (n=l, 2,..., M). Thus, the geometry 
of every initial- layer in Fig. 12 a) is determined. 

Subsequently, we mimic each spherically anisotropic initial layer by a pair of effective 
isotropic sub-layers (type-A and type-B) with equal thickness at (r„ — r„_i)/2, as shown in 
Fig. |2b). The material parameters of type-A and type-B isotropic dielectrics can be inferred 
and derived from the result for radial conductivity by Sten |f27|; 



Now, the original spherical cloak turns to be a concentric isotropic coatings. Throughout this 
paper, M = 40 and R2 = 2Ri = 2A are chosen in order that the validity of Eq. (O can be 
maintained, i.e., sub-layers have to be sufficiently thin. Alternatively one could also choose a 
smaller M associated with the corresponding parameters determined by optimization methods 
and then subdivide each anisotropic layer into many more than two isotropic sub-layers. 

When converting the «-th initial layer into its pair of effective isotropic sub-layers, we need 
to pick a specific radial position for those e, and e, on the right-hand side of Eq. Q within 
the initial layer in order to determine corresponding parameters of the two isotropic dielectrics 
(Ea, Ha) and (eg, ^b) on the left-hand side of Eq. ^ and in Fig. I^b). According to 1291 . 
such a discretization mechanism using r ~ r„ in Eq. ^ for the «-th initial layer will give a 
good compromise for both forward and backward scatterings, while retaining good invisibility 
performance. 

3. Effects of Nonlinear Transformation in Nonlinear Spherical Cloaks 

Pendry's spherical cloak lH] is a Unear one, and hence it is a straight line if one plots the trans- 
formation function r' against the physical radius r. In what follows, we introduce two classes 




(3b) 



(3a) 



of nonlinear-transformation spherical cloaks, and we discuss how to restore and improve the 
invisibility performance after discretization by choosing a proper nonlinear transformation and 
by choosing a suitable compensation scheme while discretizing the original spherical cloak into 
multilayer isotropic structures. The two types of spherical cloaks considered here are classified 
in terms of the negative (i.e., concave-down) or positive (i.e., concave-up) sign of the second 
derivative of the transformation function. It is important to reiterate that all three designs — 
linear, concave-up, and concave-down — are perfect cloaks for the exact inhomogeous design, 
and we are only considering the breakdown of invisibility when the design is discretized into 
homogeneous layers. 

3.1. Concave-Down Nonlinear Transformation 

To compress a sphere of air at the radius R2 in Q.' (original) space into a shell at the region 
Ri < r < R2 in Q. (compressed) space, we propose a class of prescribed functions 
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where x denotes the degree of the nonlinearity in the transformation. When x is very small in 
Eq. (HI, the curves are difficult to distinguish and all approach to the same limiting case when 
x^O: 

RoLnlr/Ri] 



r{x) = 



(5) 



Ln[R2/Ri] ■ 

Thus the parameters e, Jl in the transformed coordinates can be written in term of e', Jl' in 



the original space by 

e(r) =Ae'{r')A'^/det{A), J[{r) ^ AJl' {r')A^ /det{A) 
where A is the Jacobian matrix with elements defined as A^^, = dr/^/ dr'. 
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Fig. 3. The concave-down nonlinear transformation. When x < 0.1, all mapping curves 
are overlapping, meaning they effectively lead to the same performance. When x becomes 
large, the curve will nearly become a step function over IX < r <2X. 



Note that all curves belonging to those transformation functions in Eq. (HI have negative 
second derivative with respect to the physical space r, and we term this class of transformations 



the concave-down nonlinear transformation. The nonlinear transformation function in Eq. (|4]i 
only depends on the radial component r in the spherical coordinate system (r, 9, (p). Thus it is 
easy to find that the Jacobian matrix A is diagonal. Considering that the original space is filled 
with air (e' —Jl' — 1), Eq. (|6]l can be rewritten as 

e = /r = dmg[Xl Xl, Xl]/XrXeX^ = diag[^, ^] (7) 

AgAA A,-A|f, A,-A0 



where 



^r=^ = — ^TT (8a) 



:A. = ^= ^ (8b) 



denoting three principal stretches of the Jacobian matrix. Finally, the desired parameters in the 
compressed space (Ri < r < Rj) are shown to be 

R^+Ur-^-R^y 



xR] {R^2 - ^1 ; 



xR]Rl+' 



ee=^e = s,=^, = j^^-^^. (9b) 

Given such radial and transversal parameters, it is obvious that Eq. ^ cannot be simply solved 
due to the complicated anisotropy ratio Ae,,„ ~ et/£r= y?'R?^ j (r^ — )^ in our current case. In 
addition, one can also find that when the virtual space is equally discretized using a concave- 
down nonlinear transformation, the segmentation of the initial layers in Fig.|2]in the physical 
space r is 



r„(x) 



R\-R\ n 



R\ M 



■Rx, n=l, 2,...,M. (10) 



Here, we will study the effects of the nonlinear factor x in the near-field and far-field of the 
discretized nonlinear-transformation spherical cloaks. We fix M=40, and consider the cases of 
x=0.1, x=l, x=4, and x=10 in concave-down nonlinear cloaks, shown in Fig.|4] It can be seen 
from Figs. |4|a-d) that when x increases, the magnitude of electric field increases significantly 
inside the cloak. This is because more energy is guided towards the inner boundary of the 
PEC core, which in turn makes the cloaked PEC more visible to external incidences. To prove 
that the large electric fields only occur in the region Ri < r < R2 and to study the effect of 
X on individual near-field perturbations more explicitly, only the fields outside the cloak are 
presented in Figs. Ue-f). One can again confirm that, for discretized concave-down nonlinear 
spherical cloaks, smaller x leads to smaller disturbance in electric fields in the outer space. In 
particular, the degradation in the invisibility in the forward direction is proportional to the value 
of X (see Fig.Htg) and Fig.Sfh)). 

In Fig. |5] we consider the far fields of those cases of Fig. |4] For comparison purposes, we 
introduce two more curves, PEC and Pendry's cloak. It should be noted that the curve corre- 
sponding to Pendry's cloak is a straight line, and the same segmentation procedure as discussed 
above is applied to it, which is the same as equally dividing it into M initial-layers in the physi- 
cal space (see Fig. O. From Fig. |5] it can be seen that for concave-down nonlinear transforma- 
tions, the cloaking property is better retained when x is small. Except for the forward direction. 
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Fig. 4. Near-field interaction in the presence of the proposed nonlinear spherical cloak with 
different values of x for the concave-down class. The inner region < r < Ri is filled by 
PEC. M is set to be 40. Frequency is 2GHz. The total electric fields are plotted only in the 
region r > Ro in (e-h). 
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Fig. 5. Bistatic RCS of the concave-down nonlinear cloaks. M = 40. 



concave-down nonlinear transformation with small x can achieve even lower RCS over a wide 
range of observation angles. As x keeps increasing, their RCSs are increasing dramatically and 
can be larger than that of a uncloaked PEC core (e.g., x = 10). No matter how large x becomes, 
the suppressed backward scattering is still maintained |71 though it is not exactly zero due to the 
discretization, but the forward scattering will become higher and higher, which is not desired 
in the cloaking application. Those far-field characteristics are consistent with the perturbation 
in corresponding near-field patterns in Fig.|4] 

3.2. Concave-Up Nonlinear Transformation 

We can also propose another class of prescribed functions 

= (11) 

It is found that when x ^ 0, it also approaches to the same limit as Eq. (|5]l. In the same manner, 
we can obtain the desired parameters in the compressed space {R\ < r < R2) 
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Note that the case x = 1 is exactly Pendry's linear spherical cloak. 

It can be seen in Fig.|6]that all curves belonging to Eq. (fTTl i have positive second derivatives, 
and we term this class of transformations as the concave-up nonlinear transformations. Also, it 
is found that the hmiting case of Eq. (|5]l, i.e, x 0, is the dividing line between the convace- 
down and concave-up classes. Since x <C lin Eq. (fTTT i is actually very close to the situation of 
X ^ 1 for concave-down nonlinear transformations, we will not revisit that case here. Instead, 
we only consider x = 1, x = 4, and x = 10 concave-up transformations [Eq. (fT2li1. and the 
segmentation for concave-up nonlinear spherical cloaks in physical space r turns out to be 
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We are interested to see which class will give rise to better invisibility performance, after 
discretization, if all the other conditions are the same. We still plot the near field of concave- 
up cloaks at "x=l", "x=4", and "x=10". Comparing Fig. |4] with Fig. |7] one can see that: (1) 
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Fig. 6. The concave-up nonlinear transformation. Wiien x = 1, it is exactly Pendry's linear 
cloak, which can be verified by inserting x = I into Eq. l ll2t . When x gets extremely large, 
the mapping curve will nearly become a sharp impulse when r^2X. 



the peak amplitude of the electric field is proportional to x in both classes; (2) for each value 
of X, concave-up transformation cloaks have lower peak amplitudes than concave-down trans- 
formation cloaks, i.e., the perturbation inside the cloak region is smaller; (3) the perturbation 
outside the cloak of concave-up transformation cloaks is also lower than that of concave-down 
transformation cloaks in all cases of x; (4) the invisibility performance is better maintained for 
concave-up transformations even when x^ 10 or even larger. Hence, spherical cloaks based on 
concave-up nonlinear transformations exhibit better invisibility after discretization. 

Fig. [U exhibits some interesting features which distinguish it from the far fields of concave- 
down transformations in Fig. |5] Over a wide range of angles, the invisibility performance of 
spherical cloaks based on concave-up nonlinear transformations are more stable under varia- 
tions in X. More importantly, when x increases (e.g., x ~4 or x ~ 10), though the RCS will be 
a bit larger than Pendry's linear cloak at most angles, the RCS near the forward direction will 
be reduced greatly compared with that of classic linear cloak. Such far-field phenomena are 
also connected with the near-field patterns. It is important to point out that the increase of x in 
concave-up transformations [see Figs.|7];a-c)] will push the "hot" areas (where the electric field 
is very high) further and further away from the spherical PEC core, so the induced shadow and 
the far-field pattern become stable with increasing x. However, the increase of x in the concave- 
down transformations has little effect in shifting away the "hot" positions from the core [see 
Fig. Sta-d)], which results in larger interaction with the PEC core. Hence, the RCS reduction 
in the far field is degraded with the increase of x in the concave-down class as shown in Fig. |5] 

4. Justification of Improved Segmentation 

It has been pointed out that usually we require a large value of M during the discretization pro- 
cess, otherwise the conversion scheme in Fig.|2]will lose its accuracy. We also find that a nonlin- 
ear transformation with a given x of its concave-up class will yield a better invisibility by taking 
the segmentation as Eq. (T3[ in the physical space (this is also the only way to determine sub- 
layers of type-A and type-B dielectrics in practice). Here the justification is given. The concave- 



up class with x = 4, which is characterized by the mapping curve r (4) = / h {r/Ri) —I 

is considered as an example. 

In Table I, different segmentation r„ (x) sets (see Eq. [TJt for real radius r are applied to 




0.2 0,4 0.6 



Fig. 7. Re[EjQj^j] on x-z plane for the proposed concave-up nonlinear transformation cloaks 
at X = 1, X = 4 and x = 10, respectively. Figs. |7{a-c) present total electric fields in all 
areas while Figs.|7te-f) show the total fields only outside the cloak (r > R2). M = 40 and 
/ = IGHz- In order to show the disturbance in surrounding outer space, the plot range is 
larger (from -0.6m to 0.6m) in the right column. 
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Fig. 8. Bistatic RCS of the concave-up nonlinear cloaks. 



Table 1. The total scattering cross section (unit: dBsm) for concave-up nonlinearly trans- 
formed spherical cloaks (x = 4) under various segmentations in physical space. 



r„{x) in Eq. (13) 


'•«(!) 


rni2) 


rni4) 


rni6) 


r„(10) 


r'(4) in Eq. (11) 


-3.3 


-3.5 


-18.9 


-3.8 


-3.1 



. * Note that r„{l) for concave-up cloaks corresponds to the case of equally dividing the physical space. 



r'(4), e.g., the sets of jc ^ 1 (r„(l)), jc = 2 (r„(2)), jc = 6 (r„(6)) andx 10 (r„(10)). The total 
cross sections corresponding to those four sets are compared with the best segmentation set of 
its own r„(4), i.e., equally dividing the virtual space of jf = 4 into 40 initial-layers and then 
projecting onto physical space. We find that, for concave-up class x = 4, r„(4) is indeed its 
optimal segmentation among these possibihties. 

5. Conclusion 

In this paper, we discussed two novel classes of nonlinear-transformation-based spherical 
cloaks. An approximation mechanism for such classes of spherical cloaks by concentric 
isotropic sub-layers is proposed. The optimal segmentation in the virtual space is presented, 
and the role of nonlinearity in the coordinate transformation is discussed in order to provide 
better invisibility performance for the proposed classes of nonlinear spherical cloaks in both 
the near and far field. We find that concave-up nonlinear-transformation spherical cloaks exhibit 
better and more stable invisibility performance than the corresponding segmentation of linear- 
transformation or concave-down cloaks. An interesting topic for future exploration would be to 
consider a wider class of segmentation schemes, for example by using the segmentation pre- 
sented here as the starting point for a numerical optimization (e.g., Ref. Il28l for cylindrical 
cloaks) of variable segment thicknesses and material parameters. 
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